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ABSTRACT 


A method is presented by which response data can be processed 
in the frequency domain to obtain natural frequencies and modal 
damping ratios of multi-degree-of- freedom systems. Primary atten- 
tion is focused on purely mechanical systems which possess clas- 
sical normal modes. Systems that do not possess such modes, but 

T 

whose damping matrix [ C] is such that the matrix [cp] [C] [cp] , [cp] 
being the modal matrix, is not diagonal but whose off-diagonal 
elements are small, are also considered. When the power spectral 
density of the excitation is unknown, the method requires it to 
be reasonably flat. Numerical examples are given for systems 
with one and two degrees of freedom. 
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INTRODUCTION 


The problem of processing re^oonse data of an aerospace 
structure to controlled excitation or random excitation during 
flight in order to determine the dynamic characteristics of the 
structure has recently received considerable attention in the 
United States [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11] as well as 
abroad [12, 13, 14] . Prom an applications point of view, dynamic 

signature analysis is very important since it can be used to 

j 

detect flaws or to predict the onset of flutter during subcritical 
wind-tunnel or flight tests. 

In references 1, 2, 8, 10, 12 and 14 the autocorrelation 
function of the random response is utilized to obtain natural 
frequencies and modal damping ratios in the frequency domain. 

This method of modal identification assumes that the frequency 
spectrum of the random exciting forces is flat (white noise) 
over the frequency range of the modes of interest. The method 
becomes inaccurate if the excitation is not white noise and/or 
if the structure has closely-spaced vibrational modes. Also, 
errors are introduced due to statistical causes associated 
with signatures of finite duration. 

The random decrement method, introduced by Cole [1], has 
some definite advantages over the autocorrelation function and 


^Numbers in brackets designate references. 
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has been utilized by the inventor for the on-line failure detec- 
tion and damping measurements of aerospace structures. Some of 
the advantages of this method are: (1) changes of intensity of 

the input do not affect the level of the signature, which is at 
constant anplitude? (2) the signature distortion is considerably 
less, than in the case of the autocorrelation method, if the 
input spectral density is not flat. 

The work by Cole is primarily concerned with systems having 
one degree of freedom. Very recently, Chang [11] studied the 
dynamic characteristic of aeroelastic systems using randomdec 
signatures. Multi-degree-of-freedom were analyzed, but numerical 
experiments were carried out only for systems with two degrees 
of freedom (purely mechanical as well as aeroelastic) . 

The purpose of the current investigation is to develop a 
method by which response data of an aerospace structure (obtained 
during flight) can be processed in order to determine its natural 
frequencies and modal damping ratios. The analysis is to be suf- 
ficiently general to include systems other than just purely 
mechanical with classical normal modes, in order that it may be 
applied to the Shuttle Pogo Instability program. 
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ANALYSIS OF ONE- AND MULTI-DEGREE-OF- 
FREEDOM SYSTEMS WITH CLASSICAL DAMPING 


First, systems with one degree of freedom are analyzed in 
detail. Then, multi-degree~of-freedom systems possessing clas- 
sical normal modes are considered. Finally, in Section 3, purely 
mechanical, linear systems with nonclassical (non-proportional) 
damping are discussed. 


2.1 Systems with One Degree of Freedom 

Consider a one-degree-of freedom system whose motion is 
governed by the differential equation 

X + 2 DCx + Cfx - rr J f(t) (1) 

where f(t) is the random excitation, that is f(t) is a function 
representing an ergotic random process, Q is the natural frequency, 
and £ the damping ratio. The power spectral density of the (random) 
response x(t) is given by the familiar relation 

S x (uu) = H* (uu ) H (uu ) S f (a) ) , (2) 


where s x (w) is the power spectral density of the response. 


S f (<u) 
H(o>) 


is the power spectral density of the excitation. 



f 



H* (uu ) is the complex conjugate of F(w). 

Equation (2) can be rewritten in the form 

S x (uu) = | H('JU) pS f (m). 


( 3 ) 
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where 



(4) 


Note that the power spectral density (a> ) is the Fourier trans 
form of the autocorrelation R x (t) # that is 


S X M “ J R x (T)e” lujT dT # 

— 00 


where 


R (t ) 

x 


T 
» 2 

lim ~ Lx(t)x(t + i )dt 

T ~ -I 


(5) 


( 6 ) 


The autocorrelation in terms of the spectral density is given by 


V T) - £J 

•.GO 


(7) 


Equations (5) and (7) is Fourier transform pair. Analogous 
relations exist between S^(uu) and R f (t ) . 

Let X n be the finite Fourier transform of x(t), defined by 

T 

X n = IJt x(T)e " iUJnTdT ' (8) 

2 


where 



n = 1^ 2 , 3 , ... 


and T is the time interval for which the (random) response has 
been measured. The function 


S' (uu ) = X* 
x n T n n 


( 9 ) 
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is an estimate of s x (w) at the frequency uu = uu n (see, for example 

Reference 15). The spectral density S'(w ) can be obtained from 

x n 

a random response record x(t) from (9) by means of an FFT (fast 
Fourier transform) or by means of a real time analyzer. 


2.1.1 Definition and Minimization of Error Function 

It is possible to define an error function in the mean square 
sense such that minimization of this error will yield the parameters 
of the system Q and Q, Let us define the error according to 

e(n.C) = | [S (id) - s'(u, n )] 8 diD ( 10 ) 

o x x n 


Define a vector § by 


£ - t«) = {?} 


(id 


For e(Q,C) to be a minimum 


If. ■ 2 jo [S x ( w> - S >» if— 1 - 2 (12) 


3S x (w) 


Let 

Ihtn.c); = If 


(13) 


Expanding the vector {h} into Taylor series# 

l h (£ n+1 )} = l h (l n )} + ~4 f" A i n + °M A £ n H 


( A § n , ? n+1 - ? n ) 

’ rsw' 


( 14 ) 
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In order for 

l h (l n+1 )} a °> 

the increment A£ n must be obtained from 



This equation can be rewritten in the form 

{ A? n ] = [D n J [h n ], 

where 



(15) 

(16) 

(17) 

(18) 


The elements of the derivative matrix [D] can easily be 
determined from Eq. (10) (See Appendix I). Thus, Eqa (12) can 
be solved by iterating. Starting with an initial guess {? 1 }, 
[D 1 J and {h 1 } can be computed by numerical integration. Then, 
(A? 1 ] can be determined from (17) and (? 2 3 from 

{? 2 } = (5 1 ) + (A? 1 } 

The process is continued until {A? n } = {0}, i.e. until 

< s ; i = 1 , 2 
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2.1.2 White Noise Excitation 

In the above, it has been assumed that the power spectral 
density of the forcing function, S^(uu), is known. It is possible, 
however, to determine the parameters 0 and Q even if S f (uu) is not 
known, provided it is constant {white noise) or nearly constant. 

If the excitation is ideal white noise, S f (uj) = S 0 = const. 

Notice that for lightly damped systems for which | H (uj ) j 2 has a 
very steep peak in the neighborhood of Q (and is practically 
zero away from Q) , the behavior of S x (o> ) = | H(w) |°S f (uu ) depends 
almost entirely on | H (uu ) | 3 (and only slightly on S f (uu)), provided 
that S f (uu ) is flat, that is, it has no peaks in the neighborhood 
of 0. See Figure 1. 

Assuming that S^(uu) is reasonably flat, it can be set equal 
to a constant, S 0 , without appreciable loss of accuracy in deter- 
mining Q and C (for lightly-damped systems). This follows from 
the discussion of the preceeding paragraph. Equation (9) for 

the error function can be rewritten in the form 

2 

.i'Q.C) = J 0 [ S o SH(uu) p - S^(w)J dm (19) 

The actual value of is net important in the minimization of e 
and it can be set equal to unity. This is shown in Appendix II. 

2.2. Systems with N Degrees of Freedom 

The differential equations governing small oscillations of 
viscously-damped systems with N degrees of freedom are 


[mHx] + ( C] { x3 4 [KjfxJ = { F ( t ) J 


( 20 ) 
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where all symbols have their usual meaning. Let [ qpj be the 
modal matrix assumed normalized with respect to the inertia 
matrix, [M] . Then 

tcp] T [M] rep] = [I) 

and (21) 

[cp] = r» s j 


2.2.1. Systems with Classical Normal Modes 

Assume that the system, governed by (20), possesses classical 
normal modes, that is* assume that the damping matrix is such that 

[SPl T |C][?] " [R] = T2W (22) 


The identification problem for the N-degree-of-freedom system, 
analogous to the problem for the simple oscillator discussed in 
Section 2.2.1, is to determine uu^, (i = 1, 2, ...,N) from one 
or more random responses, x^(t), due to the random forcing function 
{F(t)3. 

Uncoupling the differential equations of motion (by means of 
the normal-coordinate transformation) and proceeding as in the 
case of the single-degree-of-freedom system, it is easy to show 
(see, for example. Reference 16) that 


N 


_L 

2TT 


N m 
C 


4 4 

r-1 s= 1 


. cp . 

J - 5 - | H*(«»)II ('i.)S- , (m)e 
yu * J r s ff 
r s -■* r s 


rii'i 


dm (23) 


(t)x ^ (t4T ) 
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where 


R v (t) = x. (t)x. (t+T) is the cross-correlation of the 
lj 1 J 

response 

Sf ^ (yj) = cross-spectral density of the forcing function 


i 1 


H r ( oj ) 


{£} = [^) T {Fj 

i - (i y 1 * ) 


(24) 


H* (uj ) = the complex conjugate of H r (<jj) and all other 

symbols have been previously defined. 

Since 


V- (T) = * J s x. x ,<^ )eiluTd "' 

13 -» 1 J 


(25) 


where S (<jj) is the cross-spectral density of the response, 

x i x j 

it follows that 


N N 


S (uu 

x i x j 


r* ^ir^is 

'■[ Z H;<-)H s Ms f f ( 

„ , _ - r s rs 


uu ) 


26 


r=l c=l 


or 


N N 

V 

L 


<p. Cp • S, £ (rjj » 
r ir Y ]s frf 3 ' 


S (uu) = ) ■ _ r 

X i X j , ! uu 3 -iu s -i (2C it'Ji ) I a' 3 -uu 2 + i (2C ujuu ) | 

J r=l s=l l r *r r JL s b s s'J 


Assume now that the cross-spectral density of the excitation is 
zero, i.e. 

S f f h) = 0, r ? s 
r s 


s f £ On) = s (■•-) 

r r r 


and let 
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Then, 


4 5 C ^ir Cp -ir S r (u,) 

S (u) ) = ) 3- . r , 3E , J E 

X^X . U ... 2 ...2 \2 , a + 2 ... 


(uu 2 -w 2 ) 3 +4£ 2 uj 2 w 2 
r=sl r b r r 


The finite Fourier transform of x.(t) is 


, » 2 -iQ t 

X jn = T J T X j (T)e dT ' 
2 


where 


= n (^f) * n = 1, 2, 3, . . . 


Notice that the function 


s ' „ («0 = lx. X* 

XjXj T ]n ]n 


is an approximation of S (uu ) at uu = 0 , 

X j X j n 


2.2.2 Definition and Minimization of Error Function 
Define an error function e^ according to 

,,00 

e j{uji^ 2 , •**' ^n/i/ 2, ***' ^N } = jo fS x.x. (uj)_S x.x. ^^ 2du) 


Suppose that p natural frequencies lying in the frequency band 
(f; T Q ) and the associated modal damping ratios are to be iden- 
tified by processing m records x^(t), i - 1, 2, . .., m. Following 
a procedure similar to that for a system with one degree of freedom. 


{h} . = (j = 1, 2, . .., m; r = 1, 2, ..., 2p) , (31) 

J °% r 
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Where the vector {£} is defined by 


{?} = ^1^2, 


"'p, t l, C 2. 


C p ) T (32) 


The increments {A§}j at each iteration are computed from 


where 



(33) 

(34) 


The elements of the derivative matrix in (34) can easily be derived 
by differentiation, as in the case of the simple oscillator (App. I) . 

Note that the m random records, x^(t), j = 1, 2, . . . , m, 
will gi* ^ rise to m iterative matrix equations like Eq. (33) . 

Each system of equations yields (after convergence) values for 
tu. <jj 0 . . . , uj Ci Co C • The results are then averaged. 

It should be pointed out that for light ly-damped systems 

the products | H (u) ) 1 1 H g ^uu)| are small for r / s in comparison to 

the same products for r = s. This implies that the double sum of 

Eq. (26) reduces to a single sum for lightly-damped systems even 

though the power cross-spectral densities of the excitation, 

S- - (uu), may not be zero, provided that they are reasonably 
x r r s 
flat . 
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3. SYSTEMS WITH NONCLASSICAL NORMAL MODES 

This section deals with nonclassically damped linear systems. 
Using a matrix perturbation technique, general expressions are 
derived for the correlation and spectral density matrices. Spe- 
cial forms of these matrices considering only a first-order per- 
turbation are also given. Returning to Eq. (22), let us assume 
that [R] is not diagonal, but the damping matric [C] is such that 

[R] = [D] + e[&), (35) 

where e is a small quantity. It is shown in Appendix III that 
the cross correlation matrix [R x (t)] is given by 

[R x (t)) = IT *<i#)] [S f <ou)] tM- («u)] T e iuJ ' r du U [cp] T , (36) 

— CD 

where 

(S f to)] = [cp] T [S F (u))] [tp) (37) 

The power spectral density matrix [S (uu ) ] is given by 

A 

[S to)l = [ep] [7 * (id) ) (S_ (w )] [7 (uu ) ] T [ tp ] T (38) 

A J- 

In Eqs. (36) and (37), the matrix [ Y (uu ) ] is defined by 

--1 

[ Y (uu ) ] = rH(uu)jjjI] + ieuj[fl] rH(tu)JJ , (39) 

where 

rH(uj)J = [-uu 2 [I] + iuuTDj + ruu 2 J ] (40) 
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3.1 Results for a First-Order Perturbation 

The matrix [Y {uu)] for a first-order perturbation is 

[Y(u>))i = rH(uu)J - ieuurH(uj)J [J»] rH(u))j (41) 


Substituting this expression in (36) and (38), yields 

[R x (T )Ji = [Cp] j tH*(w)J [S f (««)] rH(uu)Je atJJT dui[cp] T 


(42) 


+ [Cp] J uuH* [£)HS f - S f HjS] He 1UJT duu [Cp] T 


and 


[S x (0i)] x = [cp] rH* (uu ) j [S f (uu) ] rH(uu)j [cp] 


+ ieu) [cp] [H*SHS f - S f H* &H] [Cp] 


(43) 


Notice that the leading terms of (42) and (43) are identical to 
the cross-correlation and power spectral density matrices, respec- 
tively, of classically damped systems (for convenience, the brack- 
ets have been ommitted in the second terms of (42) and (43)). 

An error function similar to that defined by Eq. (30) can be 
set up using Eq. (38) or (43). Minimization of this function 
would yield the natural frequencies of the system, the elements 
of [D] and the elements of [£»] . For systems with nonclassical 
linear damping, however, it may be more efficient to determine 
the complex eigenvalues of the reduced system of 2N equations 
rather than usj.ng the approach discussed herein. 
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4. AN ALTERNATE METHOD OF MATCHING 

POWER SPECTRAL DEN&ITIES 

4.1 Systems with One Degree of Freedom 

A trial-and-error method can be employed to minimize the 
error function defined by either Eq. (10) or (30). The procedure, 
which is straight-forward, is described below. 

As in the iterative scheme discussed in Section 2, an initial 
estimate of Q and Q is required. The difference between the two 
methods (Newton-Raphson and trial-and-error) is that the latter 
requires initial specification only of the ranges in which Q and 
Q lie. Given these ranges, the damping ratio Q is held fixed and 
the natural frequency Q is allowed to vary within the assumed range 
until a value of Q is found which minimizes the error. The fre- 
quency is then fixed (at this value) and the damping ratio is 
varied through its range until the error is minimized further. 

The procedure is repeated until the error becomes smaller than a 
desired value. 

4.2 Multi-Deqree-of-Freedom Systems 

For a system with N degrees of freedom, initial estimates of 
a) . , £. (i = 1, 2, ..., N) or ranges within which these parameters 
lie are required. All uu . 's and C.'s except one, say uu. , are held 
fixed and id. is varied within its assumed range until the error 

K 

given by Eq. (30) is minimum. The procedure is repeated for all 
uu^ and Q. . After each sweep, a new set of ranges in uu^ and 
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(around the values determined by the previous set of iterations) 
is used and the procedure is repeated. The ranges used each time 
are made progressively smaller. 

It is clear that the trial-and-error procedure can be used 
initially and, once sufficiently good estimates of uu^ and are 
established, the Newton-Raphson iterative technique can be employed 
efficiently to yield accurate results in a few iterations. 
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5. PRELIMINARY NUMERICAL RESULTS 

Some preliminary numerical results are given in Tables 1, 2 
and 3. Table 1 shows that convergence is reached after six itera- 
tions for a one-degree-of-freedom system with parameters ft = 1.0, 

C “ 0.04. The Newton-Raphson scheme discussed in Section 2 was 
used. Two sets of iterations are shown, one with initial guesses 
ft = 0.8, Q - 0.06 and the other with ft = 1.3, C =0.06. White- 
noise excitation with power spectral density 1 was assumed. 

In Tables 2 and 3, numerical results are ,3hown for systems 
with two degrees of freedom. The parameters of the systems are: 

% = 1.0, Ci = 0.04, '\i 8 = 1.5, Cp =0.04 (Table 2); % = 1.0, 

Ci = 0.04, (JUs = 1.1, Qz a* 0.04 (Table 3). The power spectral 
densities S (uu) for these systems, assuming white-noise excita- 

x i x i 

tion with power spectral density 1, are shown in Figs. 2 and 3. 

The trial-and-error method discussed in Section 4.2 was employed 
to generate these tables. It is seen that six complete sweeps are 
sufficient for convergence. Within each sweep, each of the para- 
meters was incremented in ten equal steps (within its range) until 
the error was minimum. 
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APPENDIX I 

Elements of Derivative Matrix of Equation 17 


Consider the function H(w)| 8 in the form 


and let 


H (uu ) 


2 _ 


Q 4 


(Q 3 - uu 8 ) 8 + 4C s Q 3 uu a 


A *= (Q 2 - uu 2 ) 8 + 4C 2 rfiD 8 , 

B = ~ = 4(0 3 - Qw 3 ) + 8Quu a Q 2 , 

C = |- = 4 (317 s - uu 3 ) + 8 w 3 C®, 

D - = 8Q 2 (i> 2 C, 

E = fjj = 


The derivatives of | H (uu ) | 3 are 


To 


an 


4Cf b n * 

Vft / T - J* n 


=X * -f n3 - (x> 


^ | H (uu ) I s = - f 8 n 4 


sc 


H(-d) 


2 _ f 2D' 

* Va 3 




a" | TI , X 12 4D ^ l ZVtt £. \^4 

aoac l H{,u) l = “ F ^ + v'a 3 ' " a 3 ) n 


/ 2DB E \ 


Thus, the elements of the derivative matrix of Eq. (17) can 
be determined {by inversion) from the following derivatives of 
the error function; 
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St oo 

sfJci = 2So JoL % ^ I H ) I 


+ {% |Hto) I s - S^(ui))- a , a * ^~ |H(ui) | 3 Jdcu 


lie [Kf b .m \ s 

'iff ^JoIMa - 7? u ) 


+ (hS? „ a* t m \\.(±2tf 8BQ 3 ’C 2B 3 TV . 

+ v a s x (uj) M“a - La 2 - -t j;/*" 


d 2 e 

i?F 




+ (¥*- s x {w) )-(¥- 


d 2 e 




- 2SL f < f'ilf B n* V D -* N \ 

■ 2 So J 0 iMt * A» n A - x n i 


* (V- *>>)•(- ffl- *>)j 


do) 


Note white noise excitation of power spectral density Sq has been 
assumed, and that = Q, = r. 
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APPENDIX II 

Proof that for Ideal White Noise the Constant Value 
of the Power Spectral Density Need not be Known 


Consider a linear system with one degree of freedom excited 
by two different forcing functions f* (t) and f 3 (t), representing 
ergodic random processes. Let x x (t) and ^ (t) be the responses 
to fi and f s , respectively. If the excitations are related by 

f 2 = af x , 

where a is a constant, since the system is linear, the responses 
will be related by 


Xs - ax,, 

From the definition of power spectral density, 

S (iu ) = a 2 S (uj) 


and 


Noting that 


S f (w) = a 2 S f (uu) 
£ s r i 




¥ 


T 

2 

1 

2 


, » — 1UUT . 

Xg (t )e dT 



I 


the estimates of the power spectral densities S' (w) 

x ? 

of the responses are related by 


and S r (ui ) 

x i 


S' {'*0 = a 2 S (<») 
ah x i 
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The error functions e* and e_ , 

r i r s 

00 ft 

r 

e* = v S (uu ) - S (uu ) j duu, 

J L Xj Xi j 

o 

% = j L S *. (u,) " S Xe (,U) j duJ ' 

0 

are obviously related by e~ = a*e- . Therefore, from 

r 2 r x 


de- 

if 


i ( f?) ■ ic/) 


it is clear that the value of a does not enter in the determina- 
tion of 0 and £ (through the minimization of an error function 
defined according to equation Id, since it cancels out. 
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APPENDIX III 

Response of Linear Nonclasslcallv Damped 
Systems to Random Excitation 


The differential equations governing small oscillations of a 
viscously-damped/ linear system with N degrees of freedom, in 
matrix form, are 

[M]{x) + [C){x} + [K]{x] = {F (t) ] , (1) 


where [M] , [ C] and [K] are the inertia, damping and stiffness 
matrices, respectively, {x} is the response vector, and {F(t)) 
is a vector of applied random forces. Let (qp] be the modal matrix 
of the undamped system. If [ cp] is normalized with respect to [M] , 

[C()] T [M] [<p) = [I] 

T . (2) 

ttp] ^ [m] [cp] = rw a j 


Let 


ttP] T [C][<p] - [R] 


(3) 


The system is said to possess classical normal modes if [R] 
is diagonal. Conditions for the existence of such modes have 
been given by Caughey [16]. It is assumed herein that (Rj is not 
diagonal but that its off-diagonal elements are small. That is, 
we assume that [C] is such that [R] can be written as 

[R] * [DJ 4 €[$J, (4) 


where e is small. Introducing the normal coordinate transformation 


fx} = fcp] fpj 


(5) 
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in (1) and pre-multiplying by [sp] T , yields 

[I]fp] + CR] C p} + Cp} - [cp] T {F(t)} (6) 


Equation (6) is now solved using a matrix perturbation tech- 
nique. Assume the following expansion for the vector of normal 
coordinates : n 


tpJ 


k^O 


iP J 


(7) 


substitution of (4) and (7) in (6) leads to 

n n 

+ L fDJ + e(sl ] I elC( Pk 5 

k^O k=0 (8 ) 

n 

+ r^ B J /5 k Lp k i = {f(t)} # 

k=0 

where £ f ( t ) } = [cp] T { F ( t ) } . Equation {8) can be rewritten in 
the form 

n 

^ ({p k ) + fD]fp k ) + rw s j[p k )e k 

k=0 (9) 

n 

+ N {p k 3 * £ f ) 

k=0 


or 

[i] {ft, 3 + roj £pb } + Ipq } 

n 

+ / € k ({p k ) + f DJ £ p k ) + t"-* J Cp k J + f f Plc-i J ) * (£) 

k=l 
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Since (f) does not depend on e, 

(Pb } + tmio) + - {f } 

and n 

/ t e k ({p k ) + tDj[p k } + ^ a J{p k ] -t [A] £p k _ 1 ] ) = {0} 
k«l 


Hence, 

Ip*} + fDJ {p k } + tw 3 J fp k } - - UHp^} 


Taking the Fourier transform of (10) 


or 


L -u) a [i] + iiD (dj + r^ s jj{pb (^)J = {f(»)j 
(Pb <«*)} - fH(u))J {f (id)} , 


where 

fH(u»)J = L~uu a [ XJ + iu^DJ + fuu 2 j] 

Since 

{ft, } - iw{ft ) 

from (11) with K=l, we find 

{pi (w)] » - [HJ [C] {ft (a;)) 

= - ilt[HJ [i] ['HJ {7} 


- 1 


Similarly, 

[pte ( ) } - - tllj [&) (<w) } 

- <-iu)) s fH] [*} r K| [F.] [HJ [1 (•!!)} 
= (-i"0 a f HJ [£) f HJ j - (?(<.)] 

L J 


(10) 


( 11 ) 


( 12 ) 
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and 

{p k <(j))} = tH]j { f (<u ) J 

(13) 

Now, 

n 



[p (a)} = ^ e k {p k (uu)} 
k=0 


or 

n k 
= 1 (-iuuc ) k TH 3 j_ [^3 MU ] {f («>)} 

k=l 



C p (JU ) } = tH] ][ (-i(De) k m rH]] k {f (uu)} 
k=l 

(14) 

Considering an infinite-order perturbation (K=«>), Eq. 

(14) can 

be rewritten 

in the form 



£p} = + ieu) [*] fHJ ] *{f} 

(15) 

Substituting 

in (5) 



{ x (UJ ) } = [cp] TH(tu)J [B(uu)]“ 1 (T (uu)} 

(16) 

where 

[B(u>)] = [I] + ieuu[J9] rHo)J 

(17) 

Let 

[T (»)] = CH(u»)J [B (tu ) ] 

(18) 


Equation (15) can be rewritten as 

f p (tu ) } = ( Y (a) ) ] { f (uu ) } (19) 

Now consider the matrix 

{p(t)}[p(t + t)} T 
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Using the convolution theorem the inverse transform of fp(t)} is 

£p(t)} = j 1 [Y(X)]{F(t - X)}dX (20) 

Similarly, 


{p(t + t)} T = j {F(t + t - ?)) T [¥ (?)) T d§ 


( 21 ) 


The cross-correlation matrix of the response [R (t)J is defined by 

T P 

«2 

[R_(t)] = Xim =j T {p(t)]{p(t + t ) } T dt 
T-4- --j 

or, using (20) and (21), 

T 

[R (t)] = t’f(X)][P(t-X)}ax]| j' {F(t+T-?)} T h' (§)l T d?]at ( 22 ) 

F T->® - 00 —CO 


From (22), it can be shown that 


where 


co • 

[R (T )] = -ij 


[S f (UJ ) ] = [op] [S F (ou)3 [cpj 


£x(t)3 = (cp]{p(t)}, {x} T = [p} T [cp] T , 


(23) 


00 * 

[R x (t )3 = ^ [cp] j [¥>>)] [s f (ud ) ] [Y (uu)3 T e 1 UJT duu[cp] 


Since 
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TABLE 1 Initial Guesses and Results at Each Iteration for a 

System with One Degree of Freedom. Q “ 1.0, Q » 0.04 





































TABLE 2 Initial Guesses and Results at Each Iteration. 
= 1.0# uua =1.5# = Cs = 0.04 


Initial 

Guess 

0.8 

2.0 

0.06 

0.06 

Iteration 

No. 



Ci 

Cs 

1 

0.99200 

1.5200 

0.04560 

0.04560 

2 

0.99200 

1.5200 

0.04128 

0.04416 

3 

0.99968 

1.5008 

0.04013 

0.04186 

4 

0.99968 

1.5008 

0.04013 

0.04093 

5 

0.99968 

1.5008 

0.04004 

0.04057 

6 

1.0002 

1.4996 

0.03999 

0.04042 


Assumed Ranges: 0.56 s UU 3 , s 1.04 

1.12 £ id, < 2.60 
0.42 s Ca ,Cs s 0.78 









































TABLE 3 Initial Guesses and Results at Each Iteration. 
= 1.0, ujg =1,1, — Qa — 0.04 


Initial 

Guess 

0.8 

1.5 

0.06 

0.06 

Iteration 

No. 


U)g 

Ci 

Cs 

1 

1.0400 

1.1400 

.04560 

.04920 

2 

1.0016 

1.1040 

.03980 

.04344 

3 

1.0016 

1.1040 

.03926 

.04113 

4 

1.0016 

1.0982 

.04186 

.04045 

5 

1.0004 

1.1005 

.03991 

.04008 

6 

.99988 

1.0996 

.04005 

.04000 


Assumed Ranges: 0.056 s uuj, s 1.04 

1.05 s u> 2 s 1.95 
0.042 s Ci »Cs s 0.078 











































Lty for a two-degree 
1 — 1 . 0 , U/g ® 1 « 5 / 








